Physically adequate proper reference system of a test observer 
and relativistic description of the GAIA attitude 



Sergei A. Klioner 

Lohrmann Observatory, Dresden Technical University, 
Mommsenstr. 13, 01062 Dresden, Germany 
(Dated: February 2, 2008) 

Abstract 

A relativistic definition of the physically adequate proper reference system of a test observer 
is suggested within the framework of the PPN formalism. According to the nomenclature ac- 
cepted within the GAIA project this reference system is called Center-of-Mass Reference System 
(CoMRS). The interrelation between the suggested definition of the CoMRS and the Resolutions 
2000 on relativity of the International Astronomical Union (lAU) are elucidated. The tetrad rep- 
resentation of the CoMRS at its origin is also explicated. It is demonstrated how to use that 
tetrad representation to calculate the relation between the observed direction of a light ray and the 
corresponding coordinate direction in the Bary centric Celestial Reference System of the lAU. It is 
argued that the kinematically non-rotating CoMRS is the natural choice of the reference system 
where the attitude of the observer (e.g. of the GAIA satellite) should be modeled. The relativistic 
equations of rotational motion of a satellite relative to its CoMRS are briefly discussed. A simple 
algorithm for the attitude description of the satellite is proposed. 
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I. INTRODUCTION 



Future space astrometry projects like GAIA [1-3] and SIM [4] are expected to attain 
an accuracy of 1 microarcsecond (yuas) for positions of remote celestial sources. This high 
accuracy requires general relativity to be used for data modeling. A relativistic model of 
positional observation with microarcsecond accuracy involves many subtle details. Recently 
a number of such models have been suggested (see [5-8] and references therein). The purpose 
of this paper is, first of all, to provide a relativistic definition of the physically adequate local 
(proper) reference system of a test observer (e.g., a satellite). According to the nomenclature 
adopted within the GAIA project [9] this reference system is called Centcr-of-Mass Reference 
System (CoMRS) below. As a co-product, the calculation of the observed light direction 
as adopted in [7] is explained in detail and explicitly justified. The CoMRS is intended 
to be used to describe physical phenomena located within the immediate vicinity of the 
observer (e.g., the rotational motion of the satellite, the process of observation, etc.). This 
reference system can be used to define the attitude parameters of the satellite. In order to 
define the CoMRS we make use of the parametrized post-Newtonian (PPN) version of the 
relativistic framework adopted recently by the International Astronomical Union (lAU) for 
the use for high- precision astrometry, celestial mechanics, geodesy and metrology [10]. The 
lAU Resolution B1.3 (2000) adopted at the XXIV General Assembly [10-12] of the lAU 
specifies a global reference system, the Barycentric Celestial Reference System (BCRS), and 
a physically adequate local geocentric reference system, the Geocentric Celestial Reference 
System (GCRS), in the framework of the post- Newtonian approximation of general relativity. 
Below it is argued that a simplified version of the GCRS constructed for a massless observer 
as a central body can be used as a physically adequate CoMRS. 

The problem of constructing a physically adequate proper reference system of a massive 
body (e.g. Earth) in the first post-Newtonian approximation has been thoroughly discussed 
by several authors. In the framework of general relativity two advanced formalisms are 
available. One formalism is due to Brumberg and Kopeikin [13-17] and another one is due 
to Damour, Soffcl and Xu [18-21]. For the gravitational A^-body problem both formalisms 
introduce a total of + 1 different coordinate systems: one set of global coordinates {t, a;*) 
and one set of local comoving coordinates (T, X") for each body. Note that in this context a 
body is just a material subsystem at the boTindaries of which the energy-momentum tensor 
vanishes. In the local coordinates the metric tensor possesses the following two properties: 

A. The gravitational field of external bodies is represented only in the form of a relativistic 
tidal potential which is at least of second order in the local spatial coordinates and 
coincides with the usual Newtonian tidal potential in the Newtonian limit. 

B. The internal gravitational field of the central body coincides with the gravitational 
field of a corresponding isolated source provided that the tidal infiuence of the external 
matter is neglected. 

These two requirements can be simultaneously satisfied in general relativity as has been 
shown in the framework of the Brumberg-Kopeikin and Damour-Soffel-Xu formalisms. It is 
clear that this fact is closely related to the validity of the Strong Equivalence Principle in 
general relativity. These two formalisms complement each other by elaborating the theory 
from slightly different points of view. The formalisms deliver (1) an elegant description of 
metric tensors in both the global and local coordinates and the closed-form transformations 
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between them, (2) an improved description of the structure of the gravitational field of 
each body by means of a set of its multipole moments which are linked in an operational 
way to what can be observed in the local gravitational environment of the body, (3) a 
description of the influence of the external gravitational flelds in the local reference system 
by means of some suitably defined tidal moments, (4) post-Newtonian translational and 
rotational equations of motion of the bodies with full multipole structure, (5) physically 
adequate equations of motion of a test particle in a local reference system, (6) physically 
adequate relativistic models for many kinds of observations (VLBI, high-accuracy positional 
observations, remote clock synchronization, etc.). The lAU 2000 Resolution B1.3 [10] that 
defines the metric tensors of both the global BCRS and the local GCRS, is based on the 
results on these two approaches. 

In the framework of the PPN formalism with two parameters /? and 7 the theory of phys- 
ically adequate local reference systems was developed in [22, 23]. It is clear that because of 
possible violation of the Strong Equivalence Principle in some alternative theories of gravity 
it is, generally speaking, impossible to construct a local reference system possessing both 
properties A and B. In [23] it has ben shown how to construct local reference systems which 
possess either property A or property B. It has been also demonstrated that for relativistic 
modeling of astronomical observations one should normally prefer the local reference system 
with property A. The theory of local PPN reference systems with the PPN parameters has 
been developed in [22, 23] with the same degree of details as it was done in general relativity. 
For the limit of general relativity /9 = 7 = 1 the formulas of [23] coincide with those of the 
Brumberg-Kopeikin and D amour- Soffel-Xu formalisms as well as with the formulas from the 
lAU 2000 Resolutions. 

The problem of defining physically adequate local coordinates for a massless body (test 
observer) is much more simple than the problem for a massive body (e.g. for Earth). This 
problem has been discussed many times in the literature. Let us mention, for example, 
the work of Ni and Zimmermann [24] where a local reference system of an accelerated 
observer has been constructed explicitly up to the terms of third order relative to the local 
spatial coordinates. As it has been noted in [17], the results of the Brumberg-Kopeikin and 
Damour-Soffel-Xu formalisms can be directly applied to define in an elegant way a physically 
adequate local reference system of a massless body. Indeed, one should just consider the 
limit where the gravitational potential of the central body vanishes. It is clear that the same 
procedure can be applied also to the local PPN reference system [23]. Exactly this will be 
done below. The resulting reference system represents a natural choice for the physically 
adequate CoMRS of a massless observer. This reference system can be applied to model 
physical phenomena in the immediate vicinity of the observer. Two examples will be given 
below: the relation between the observed direction toward a light source and the relevant 
BCRS parameters of the light ray, and the rotational motion of the test observer (satellite). 

Let us summarize the most important notations used throughout the paper: 

• G is the Newtonian constant of gravitation; 

• c is the velocity of light; 

• 13 and 7 are the parameters of the parametrized post-Newtonian (PPN) formahsm 
which characterize possible deviation of the physical reality from general relativity 
theory (/3 = 7 = 1 in general relativity); 

• the lower case latin indices i, j, k, .. .take values 1, 2, 3; 
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• the lower case Greek indices a, (3, .. .take values 0, 1, 2, 3; 

• repeated indices imply the Einsteinian summation irrespective of their positions (e.g. 
a' U = 6^ + a? 6^ + b^, a'^ ba = a° 6o + bi + a^b2 + 63); 

• parentheses surrounding a group of indices denote symmetrization, e.g., ^j(jfc) = 

• brackets surrounding two indices denote antisymmetrization, e.g., A^ykj — 

• a comma before an index designates the partial derivative with respect to the corre- 
sponding coordinates: A = dA/dx'^, A^i = dA{t,x)/dx^; for partial derivatives with 
respect to the coordinate times t and T we use special notations = dA{t,x)/dt 
and A^T = dA/dT- 

• a dot over any quantity designates the total derivative with respect to the coordinate 

dA 

time of the corresponding reference system: e.g. A — ——. 

at 

Sections II, III and IV are devoted to the definitions of the metric tensors of the BCRS 

and the CoMRS, and the transformations between these two reference system, respectively. 
The tetrad induced by the CoMRS coordinates at the CoMRS origin is discussed in Sec- 
tion V. Section VI elucidates the equivalence of several possible ways to calculate the 
observable direction toward a light source from the relevant coordinate quantities defined in 
the global BCRS. The post-Newtonian equations of rotational motion of a satellite relative 
to the CoMRS are briefly discussed in Section VII. In Section VIII it is argued that the 
kinematically non-rotating CoMRS represents a natural choice of a reference system where 
the attitude of the observer (e.g. of the GAIA satellite) is modeled. A summary of the main 
results are given in Section IX. 



II. THE PPN METRIC TENSOR IN THE BCRS 

Let us consider an isolated system of gravitating bodies. It is clear that the space-time 
is asymptotically fiat and can be covered with a single global coordinate system = {ct, a;*) 
where 

lim Qi^v^Vnv, (1) 

t— const 

g^y being the metric tensor in the global coordinate system. Here r]^,, = diag(— 1, -|-1, -|-1, -|-1) 
is the Minkowski metric tensor. In the framework of the PPN formalism [25] with two 
parameters (3 and 7 the metric tensor Qi^^, in the global reference system can be written as 
[22, 23] 



2 . . 2 



^00 = -l + ^w{t,x)-^pw\t,x) + 0{c-'^), 



c 



2(1 + 7) 

^3 



9oi = - "'^ ' " w\t,x) + 0{c ^), 
2 



= 5,,- (1 + -7^(t,a;)) +C»(c-^), (2) 
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where Sij is the Kronecker symbol. Here and below 3-dimensional coordinate quantities 
( "3-vectors" ) referred to the spatial axes of the corresponding reference system are set in 
boldface: e.g. x = x^. A harmonic-like gauge for the global PPN metric tensor is adopted 
here. Precisely speaking, the global metric tensor satisfies the usual harmonic gauge {g = 

% {{-gfV') - (3) 



dx' 

in case of general relativity (3 — ^ — 1. This requires 

w,t + wl^ = 0{c-^). (4) 

In accordance with the standard PPN framework as described, e.g., in [25] the metric po- 
tentials w and are assumed to obey the equations 

w,u - ^ w,tt = -4 TT G a + 0{c-^) , (5) 
w'jj^-AnGa' + Oic-^), (6) 

where 



a = 1 (^T°° + 7T'='= + iT°°(37-2/3-l)w) +0(c-^), (7) 

a' = ^T°^ + C»(c-2). (8) 

Here, T^" arc the components of the energy- momentum tensor in the global reference system 
and w in (7) is needed only to Newtonian order where it coincides with the Newtonian 
potential. Because of requirement (1) the solution of (5)-(6) can be written in the form 

w^'it, x)^gJ ^^^^^V + ^G^ I a^{t, x')\x - x'\d'x' + Oic-'), (9) 

where = w and = a. It is clear that this formulas for w and together with the 
Newtonian continuity equation 

a, + a], = 0{c-') (10) 

are compatible with the gauge condition (4). The metric (2)-(9) is equivalent to the PPN 
metric with coordinates [t^^.x]^^) in the standard post-Newtonian gauge as used, e.g., in 
[25] up to a trivial gauge transformation 



where x is the superpotential 

X = ^G J a{t, x') \x - x'\ (fx' + 0{c-^), (12) 
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so that 



X,ii = w + 0{c-^). (13) 

For /9 = 7 = 1 the definition (2) and (7)-(9) coincides with the definition of the BCRS 
[10, 11] given within general relativity. 



III. THE PPN METRIC TENSOR IN THE COMRS 



The CoMRS is a physically adequate reference system of an observer (T, A"") the mass 
of which is so small that its infiuence on the background space-time can be neglected. The 
metric tensor in the CoMRS can be derived from the metric tensor of the GCRS by setting 
the gravitational potential of the central body to zero. Below we modify in this way the 
PPN version of the GCRS as constructed in [23] . Again by setting ^ — P — lin the formulas 
below one can restore the formulas which could be derived directly from the GCRS adopted 
by the lAU [10-12]. The metric tensor in the CoMRS reads 

Goo = -l + ^W{T,X)-^(3W'(T,X) + 0(c-'), 

Goa = -^^^^W%T,X) + Oic-'), 

Gab = 5a, (l + ^7 >V(r, AT)) + 0(c-^), (14) 

where 



W{T,X) = Q„(r)A'" + WT(r,A'), (15) 
>V»(T,Ar) = Ua,cC,{r)X- + W!i{T,X), (16) 

where eijk = {i — k){k — i)/2 is the fully antisymmetric Levi-Civita symbol. Potentials 

Wt(T, X) and yV^{T, X) are the external tidal gravitational potentials (both arc (9(|A'p)) 
which describe the manifestation of the external gravitational field in the CoMRS. The 
terms Qa{T) X"' and \eabcCb{T) X^, Qa and Ca being arbitrary functions of time T, are 
hnear relative to \X\ and describe the translational and rotation motion of the CoMRS. 
The Qa is the acceleration of the momentarily co-moving locally inertial reference system 
relative to the CoMRS origin. In other words, an accelerometer placed at the CoMRS origin 
measures —Qa (see. Section VIII of [23] where the equations of motion of a test particle 
relative to the local PPN coordinates were derived). If the observer (satellite) is a drag- 
free satellite one can set Qa = 0. If the observer (satellite) is equipped with some kind of 
thrusters, Qa{T) can be used to describe their infiuence. Non-gravitational forces can be 
also described by choosing some special model for Qa(T). In the following we consider Qa 
as arbitrary function. The Ca defines the rotational motion of the spatial axes of the CoMRS 
relative to the momentarily co-moving Fermi- Walker transported locally inertial reference 
system. Clearly, the equations of test particles relative to the CoMRS with Ca 7^ contain 
Coriolis forces. Possible choices of Ca and its relation to the rotational matrix 7?."^ in the 
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transformation between the CoMRS spatial coordinates and the spatial coordinates of the 
BCRS will be discussed below. 

Here the harmonic gauge conditions are again assumed to be valid in the /3 = 7 = 1 limit. 
This imphes an equation similar to (4) for potentials W and W" that in turn gives 

Q„ X- + >VT,r + >V^,„ = 0{c-^). (17) 
Now from the results of [22, 23] one gets 



yVT(T,A')= w{t,x)-w{xo)-w^j{xo)ri 

+i (-2(1 + 7) vl {w% X) - w\x,) - wy{x,) ri) 

+ (1 + 7) ^^o - w{Xo) - w,j{x^) ri) 

+ (1 + 7) w''^{Xo) rlri + ^jw{Xo) + (^ - /? - 7) Krlf 
+ (1 - 2/3 - 27) Q„ nt rl ai - 7 < w, {x,) ri 
+ lirlQan\al + ^{-f-2)rldlrl 
+ 2(1 - /3) [w{Xo) + a>;] (w{t,x) - w{xo) - wj{xo)rij 



+0{c- 



(18) 



W^(r, X) = n% \w% x) - w\x,) - w%x,) ri - vi {w{t, x) - w{Xo) - Wj{Xo)ri) 

Here rl — — xl{t), xl{t) are the coordinates of the origin of the local reference system rel- 
ative to the global one, and = dx^^/dt and a\ = d'^x'^^/dt^ are its velocity and acceleration, 
respectively. For any function of A{t, x) we use the shorthand notation A[Xo) = A{t, Xoit)). 



IV. TRANSFORMATION FROM THE BCRS TO THE COMRS 

The coordinate transformations between the BCRS and the CoMRS read [22, 23] 

^ = ^ - ^ (-4 + < <) + 1 + + rl ri + C{t, x)) + 0(c-^), (20) 
A"" = 7^^. [ri + i {^-vi v'; + V^^ r^ + V^^' r^ r^)) + 0(0"^), (21) 

A{t) = ]^vl^w{x,), (22) 
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m = -I vt + 2(7 + 1) vlw\x,) -{l+D vlw{x,) + (/? - ^) ^'(^o), (23) 

--\vlvi + 2(1 + 7) w% x^) - (27 + 1) vi ^(Xo), (24) 

B'^{t) = -v^: 7^«) Q" + (1 + -f)w^''^\x,) - ^vi' w^^\xo) + ^jS'^ w{xo), (25) 

C{t,x) = ^{^-2)rl{airi), (26) 

V'^{t)^S'^^w{xo), (27) 

V'^\t) = ^ 7 ((5^^a^ + (^^'^a^^ - S^'^al) . (28) 

Besides that, the rotational matrix TZ^ from the transformation of the spatial coordinates 
(21) is related to Ca from (16) as 



c^n^TZ'^^ = -(l + 7)%•fc7^\c, 

+(1 + 27) vlwj]{xo) - 2(1 + 7) w^''^\xo) + ^;i*7^^] Qa + 0{c-'). (29) 

The dynamically non-rotating version of the CoMRS characterizes by = (i.e. no Coriolis 
forces in the equations of motion of test particles) and, as it follows from (29), has specific 
spatial rotation relative to the BCRS consisting of geodetic (de Sitter), Lense-Thirring and 
Thomas precessions. The GCRS of the lAU is defined to be kinematically non-rotating, 
that is, it has no rotation of spatial axes relative to the BCRS (i.e. 7^° = 5" in this 
case). Although Coriolis forces appear in the equations of motion relative to a kinemati- 
cally non-rotating reference system, this choice is especially advantageous for modeling of 
astronomical observations, since no additional orientation-related re-calculations (e.g., of 
planetary ephemeris data) are necessary. The same arguments apply to the CoMRS: the 
kinematically non-rotating CoMRS is the most convenient choice of the orientation of the 
local coordinates. For the kinematically non-rotating CoMRS TZ^ = 5f and Ca has some 
specific non-zero value defined by (29). Below we retain TZf in the formulas, but it should 
be chosen to be identity matrix 6f. 

Matching of the CoMRS and BCRS metric tensors allows one to derive also the equations 
of motion of the CoMRS origin as well. The BCRS acceleration of the CoMRS origin (that 
is, the acceleration of the observer's center of mass) reads 

< = w,i{xo)+Aai 

+ ^ (2(1 + 7) w'{xo) + (74 - 2(7 + /3) w{xo)) w,i{xo) - (27 + 1) vi w{xo) 

- 2(1 + 7) viw^^,{Xo) - viviwjixo)^ + 0{c-'), (30) 

where 
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A< = -n'; Qa {S^^ - ^ (vl 5'^ + (2 + 7) w{x,) 5^^ + \vl v^)j^ (31) 

is the BCRS coordinate acceleration of the CoMRS origin relative to the co-moving locally 
inertial reference system. Clearly, Aa^ is proportional to Qa and comes just from the re- 
calculation of the CoMRS-defined Qa into the BCRS (see, the discussion of Eq. (4.29) of 
[23]). If the acceleration AaJ, is neglected, Eq. (30) coincides with the equation of time-hke 
geodetic in the metric (2). If the gravitational fields of all N bodies can be described only by 
their masses (no further multipole moments of the gravitational field in the corresponding 
local reference system of each body), Eq. (30) produces Eq. (3) of [7]. 

From (20) with (22) and (23) follows that the CoMRS coordinate time T at the CoMRS 
origin A"" = (this is equivalent to = 0) coincides with the proper time r of the test 
observer. 



V. COORDINATE-INDUCED TETRAD FOR THE ORIGIN OF THE COMRS 



Let us construct a tetrad (e.g., [26]) co-moving with the observer. Let us first introduce 
four vectors e^^^ attached to a point on the worldhne of the observer. Here index a is the 
tetrad index which runs from to 3 and numerates the vectors. Index // is a normal tensor 
index which can be lowered and raised by contracting with the metric tensor 



6(a) li. — Qixv ©(a) , 

efa) = r^^a).. (32) 
The four vectors are required to have the property 

e^ec/?)/^ = (33) 

This equation implies that the vectors are orthogonal to each other, that vector e^Q^ is unit 
and time-hke, and that e^-^ are unit and space-like. Four additional vectors e^"^ (with tetrad 
index written as superscript) are then defined by 

e(")^e(^), = <5^, (34) 

5^ = diag(l, 1, 1, 1) is the 4-dimensional Kronecker symbol. From (33) and (34) it is easy 
to show that the conversion between vectors e^^^ and e*^'*^'^ can be performed simply by 
contracting with the Minkowski metric 



= r/^^^ef^), (35) 



where r/"^ = diag(— 1, 1, 1, 1) is the matrix inverse to rjap. With the help of these vectors 
one can represent the metric tensor at the considered point of space-time as 

5^. = e(r)e(a). = 77a;3e(r)ef , (36) 
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so that 



ds^ = 7]a/3 dx^^'Ux^^l (37) 

with 

^^(a) ^ g(a) ^3g) 

Eq. (37) shows that dx^°'^ can be interpreted as observable infinitely small time intervals 
and distances in the infinitesimal neighborhood of the considered space-time point. For a 
tetrad co-moving with the observer the unit time-like vector e|^g^ can be chosen to coincide 

with the 4- velocity of the observer, so that the projection dx^^^ = e^^^ dx^ coincides with the 
interval of the proper time of the observer dx'^'^^ = dr between the events with coordinates 
and x^ + dx^ both lying on the worldline of the observer. This means 

(39) 

The vectors e*^*) are then constrained by (36) up to arbitrary spatial rotation. This means 
that if is a solution of (36) then 

ef = ^'e(f) (40) 

with arbitrary orthogonal matrix is also a solution. 

Considering the BCRS metric tensor (2) one gets from (39) 

+ ^ (^° + (^ + ^)^o^(^o) + (/3-^)«;^(a;o))+0(c-^), (41) 
ef = [-\^l^l-i2^ + l)w{x,)vl + 2{l + ^)w\xS)+0{c-'). (42) 

The following partial solution for elf-* is then possible: 

ej"^ = -\vl (l-r^,{^-vl-riw{xS))-rO{c-% (43) 
ej") ^ 5"' + ^ (]- < < + 7 w^(aJo) 5"') + 0{c-^). (44) 



This solution for e-"^ shows that we have chosen the space-hke vectors of the tetrad is 
such a way that they show no spatial rotation relative to the spatial axes of the BCRS. It 
is easy to see from (20)-(28) that this tetrad is the coordinate basis of the kinematically 
non-rotating CoMRS (i.e. of the CoMRS with TZ^ = 6^) on the worldline of its origin 



eir^ = ^ ■ (45) 
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Therefore, the CoMRS implies also adopting a particular tetrad co-moving with the observer. 
Tetrad (41)-(44) is induced by the CoMRS coordinates at the origin of the CoMRS in the 
sense of Section 3.4.2 of [27]. This tetrad can be used to model certain kind of observables 
(e.g., proper directions as directions relative to the tetrad (41)-(44)). However, the CoMRS 
is more than just a tetrad and allows one to use all the power of the theory of local reference 
systems as mentioned in Section I. 

If one adopts a dynamically non-rotating CoMRS with = and TZ°^ defined by (29), 
the corresponding tetrad (45) will be Fermi- Walker transported, so that the Fermi rotation 
coefficients of that tetrad vanish. This is, however, an unnecessary complication for space 
astrometry, where tetrad (41)-(44) and the kinematically non-rotating CoMRS are more 
convenient. 

The tetrad (41)-(44) is written for arbitrary velocity of the observer v^. The tetrads 
e^") and e^") defined by (41)-(44) with two different velocities vl and v^, respectively, are 
related to each other by a Lorentz transformation plus additional spatial rotation V^j of 
space-like vectors: 



= (46) 
e» = P^A^ei"), (47) 

where 



K - r, 

A^ = --u'T, 
c 

/ 1 \ -1/2 

r = (^l--zy^zy'=j , (48) 

and 

n- = S'^ + ^ {K vi - vi vl) + 0{c-'). (49) 

Here, the parameter of the Lorentz transformation is the BCRS velocity of the second 
observer as seen by the first observer having the BCRS velocity vl which can be calculated 
as 



c 



el*) dx^ 



,(0) 



dxt" 



%-%+ 2 



^0 +'^k'v^/c 

^ ^ 'vi i< {vi vi)vi+l {vi vi) vi 



+ {l + ^)w{xM-vl)]+0{c-'). 



(50) 
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In the limit of special relativity Eq. (50) coincides with the special-relativistic velocity 
composition law. On the other hand, if the BCRS velocity of the first observer vanishes 
{vi = 0), Eq. (50) reproduces Eq. (12) of [7]. 

The reason for the appearance of the additional spatial rotation is the well-known fact 
that two subsequent Lorcntz transformations without spatial rotation are equivalent to the 
one Lorentz transformation with spatial rotation. This additional spatial rotation and its 
consequences for kinematically non-rotating astronomical reference systems are discussed 
in [28, 29]. The matrix (49) is equal to identity matrix 6'j if vl = 0, so that ii vl — 
the transformation between e^^^ and e^^^ is a pure Lorentz transformation with parameter 

— vl {1 + c""^ (1 -I- 7) w{Xo))+0{c~'^). Using the Lorentz transformation in its closed form 
allows one to get the expressions for the tetrad in the first post-Minkowskian approximation, 
that is, retaining terms of first order in G and of all orders in \vl\/c. 



VI. OBSERVED DIRECTION OF THE LIGHT PROPAGATION 

Let us compute explicitly the relation between the unit coordinate direction of light 
propagation n* in the BCRS and the observed direction to the light source relative to the 
tetrad (41)-(44). Let x'^{t) be the coordinates of the photon parametrized by the coordinate 
time t. It is clear that the observed direction should be defined with respect to the tetrad 
vectors ef"^ as 

where = ^ ^ is the coordinate hght velocity at the point of observation x'^ — x^. The 
differentials of x^ are calculated here along the light ray x^{t) at the point of observation. 
However, because of (45) the direction s*^"^ coincides with the CoMRS velocity of the light 
propagation at the origin of the CoMRS (that is with the velocity transformed into the 
CoMRS using the coordinate transformation (20)-(21)): 



Here Xp{T) is the worldline x^{t) of the hght ray expressed in the CoMRS coordinates. 
All partial derivatives as well as the differentials dX^ are calculated at the point of 
observation. Let us also note that (51) and (52) have one more interpretation. Let us 
consider two 4- vectors a'^ and b'^ and an observer with 4- velocity = Let us also 

car 

assume that both vectors a^^ and If' are not equal to Au^^, A being a constant. It is well 
known (e.g., [25, 27]) that by projecting each of these vectors into the observer's rest space 

and calculating the normalized scalar product of the projected vectors with respect to the 
metric g^i, one gets the cosine of the angle 9 between these two vectors as measured by the 
observer: 



a" 



P 0" 



(53) 
(54) 
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C0s9a=-jQy^. (56) 



COS 6' = TTTT, (55) 

where Pf^^, = g^^ + is the projection operator into the satelhte's rest space. Using 

(39) and (36) one can write P^j, = g^^^ + e^°^ e^^-* = 5ab e^^-*. Therefore, the cosine of the 

observed angle 6a between the incoming hght ray with the wave vector k^^ = = (l,p*) 
and a space-hke vector of the triad e^"^ can be calculated as 

Note that P'^e*^")^ = e*^")° and c^^^^e^^^^ = 5"^ and, therefore, vectors e^^") already lie in the 
observer's rest space and are normalized to unity. In (56) we used also that according to 
(33) and (39) MaC*^"^" = for any a = 1, 2, 3. This technique to compute the cosines of the 
observed angles has been used in a slightly different form e.g. in [5, 6]. It is, however, clear 
that this technique is equivalent to the two above-mentioned ways to derive s^"^ and the 
components of s*^"^ can be easily related to cos 6*0. Indeed, one has 

The difference in the sign between cos 9a and s*-"-* reflects the fact that s'-"^ is the direction 
toward the source while cos 9a characterizes the opposite direction of light propagation. 

Now one can substitute (41)-(44) into (51) or (52) and expand the denominator into 
powers of to get the explicit relation between s*^"-* and p\ The absolute value of the 
coordinate light velocity can be calculated from the fact that the light follows a null geodetic 
which means that in a reference system with metric tensor ga/s vector satisfies the equation 

9 1., k'' = ^00 + 2^oi p' + p'p^^O. (58) 
Substituting the BCRS metric (2) into (58) one gets 

IpI = 1 - ^ (1 + 7) w{x,) + ^ 2 (1 + 7) w'{x,) + C»(c-^), (59) 

where |p| = {p^ ^-p^p^ -t-p^p^)-*^/^ is the Euclidean norm of p\ Combining (51)-(52) with 
(41)-(44) and (59) one gets 

s — —n + - n X iXo X n) 

c 

+ ^^{n ■ ±0) n X X n) + ]^Xo X {n X x^) | 

+ ^ I (("■ • ^of + (1 + 7) ^(cco)) n X (cco X n) + ^ (n • x^) x^ 
+ 0(c-^), (60) 
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Here s = s"", n = n'^, and for any a* and 6* the Euclidean scalar and vector products are 
denoted as a ■ b = dijd^V = a^U and (a x b)* = Sijka^h^^ respectively. Eq. (60) coincides 
with Eq. (7) of [7]. The discussion of the Lorentz transformations of the tetrads e^^^ at the 
end of Section V allows one to conclude that Eq. (60) can be re-written as a closed-form 
Lorentz transformation (see Section 5 of [7] for further details). 

VII. RELATIVISTIC MODELING OF THE ROTATIONAL MOTION OF THE 
SATELLITE IN THE COMRS 

In principle, one can consider the rotational motion of the satellite relative to the CoMRS. 
To this end, the post-Newtonian equations of rotational motion of a satellite relative to the 
CoMRS are necessary. These equations have been derived in [20] in the framework of general 
relativity and then extended to the PPN formalism in [22, 23]. The final multipole-expanded 
form of the equations is given in Section IX. F of [23]. These are differential equations for 
the post-Newtonian spin (total angular momentum) 5" of the satellite: 

S" = (61) 

where C"" is the post-Newtonian external torque, that can be computed from the mechanical 
properties of the satellite (inertial moments, etc.) and the ephemeris data of the Solar 
system bodies (their BCRS positions, velocities, etc). 

For a satellite on a heliocentric orbit with the semi-major axes close to that of the Earth 
orbit, the largest relativistic effect in its rotational motion relative to the kinematically 
non- rotating CoMRS is clearly the geodetic precession which is of the order of ~ 2"/cty 
2 //as/hr. 



VIII. ATTITUDE DESCRIPTION OF THE GAIA SATELLITE 

It is, however, clear that these dynamical equations of rotational motion (at least the tiny 
relativistic corrections) play no role in the accurate attitude determination of the satellite. 
As it is discussed, e.g., in [2] the attitude of the satellite will be determined together with 
the astrometric parameters of the sources from aposteriori processing of the observational 
data. The attitude parameters are the parameters of the rotational matrix P"^ relating the 
CoMRS spatial axes X"" to the spatial axes X"" of the reference system in which the satellite's 
body is fixed (the latter reference system is called Scanning Reference System (SRS) in [9]): 

= P% XK (62) 

The directly observable quantities for the missions like GAIA (i.e. for scanning astrometric 
satellites) are the coordinates of the sources in the SRS tagged with the corresponding 
time of observation in the satellite's proper time 

-a^pa^b^ (63) 

where is the vector defined by (51)-(52). The observables s° should be first transformed 
from the SRS into the CoMRS with the aid of matrix and then the relativistic model 
as described e.g. in [7] should be applied to get the catalog positions of the sources. The 
matrix is clearly time-dependent and should be determined from the same observations 
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(a rough estimate of the matrix is provided by an apriori dynamical modehng of the sateUite 
rotation). The matrix P"^ can be parametrized with some Euler-hke angles or in any other 
suitable way. 

In principle, any orientation of the CoMRS spatial axes can be used to accomplish this 
data processing scheme and thus determine both the astrometric parameters of the sources 
in the BCRS and the orientation of the satellites's body relative to the CoMRS. However, the 
kinematically non- rotating CoMRS represents a natural choice of the orientation of the local 
coordinates. Indeed, in this case the difference between the CoMRS positions of the sources 
and the catalogue positions comes from a number of well-understood effects like proper 
motion, parallax, light deflection (all calculated in the BCRS) and aberration (calculated as 
discussed in Section VI above). One can also argue that Eq. (60) takes its simplest form for 
the kinematically non-rotating CoMRS: for any other spatial orientation of the CoMRS the 
formula relating the observed direction to the source s to the BCRS direction n differs from 
(60) by additional spatial rotation which exactly vanishes for the kinematically non- rotating 
CoMRS. 

One could in principle construct the tetrad directly for the SRS 



ejr^ = ne?\ (64) 

where e^") is the tetrad of the CoMRS defined by (41)-(44). This kind of tetrads was 
discussed in [5, 6]. The tetrad cj""-' can be used to compute the observed light direction 
in the same way as described in Section VI. This way is totally equivalent to using first 
the CoMRS tetrad to compute and then converting into by using (63). It is clear 
from (64) and, e.g., (51) since the contraction is associative and in particular Ple'^^^ dx^ — 
{P^e^^^) dx^ — PI {e'^^^ dx^) . In our opinion, however, keeping the matrix P^ in single 
formula (63) allows one to separate and simplify the relativistic and attitude models. 



IX. CONCLUSIONS 

The lAU 2000 relativity framework [10-12] provides not only a reasonable barycentric 
reference system for the whole solar system (BCRS) and a physically adequate geocentric 
reference system for the Earth neighborhood (GCRS). The lAU framework provides also 
suitable tools to model any kind of astronomical observations. As it was explicitly demon- 
strated above the same technique as used to construct the GCRS can be applied to define a 
physically adequate proper reference system of a test observer called CoMRS above. That 
reference system can be used to model observation processes of any kind. The coordinate 
basis of the CoMRS at its origin coincides with a particular form of tetrad comoving with 
the observer. This means that the CoMRS description of observables coincides with the clas- 
sical tetrad representation in cases where a tetrad is sufficient for modeling. On the other 
hand, the CoMRS, being a well-defined 4-dimensional reference system, is much more that 
just a tetrad and can be used to model physical phenomena where spatial extension of the 
observer plays a role (e.g., its rotational motion). In principle, any spatial orientation of the 
CoMRS is possible. However, to model astronomical observations it is preferable to adopt 
the kinematically non-rotating CoMRS having no spatial rotation relative to the BCRS. 
This choice imphes simplest possible models for observables. The attitude of the observer 
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can be then described by a 3-dimensional spatial rotation in the CoMRS. This allows one to 
separate the relativistic model and the attitude model and to simply both of them as much 
as possible. 
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